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We present a status report of our ongoing effort toward a precision determination of the kaon B parameter 
with the Kogut-Susskind quark action in quenched QCD. Results for Bk so far accumulated at /3 = 5.85, 5.93, 
6.0 and 6.2 corresponding to = 1.3 — 2.6 GeV do not exhibit the theoretically expected O(a^) behavior, but 
are apparently more consistent with an 0{a) behavior. 



1. Introduction 

The kaon B parameter is one of the crucial 
hadron weak matrix elements for constraining the 
CP violation phase of the Cabibbo-Kobayashi- 
Maskawa matrix from experiments. A pioneering 
study was carried out by Gupta, Kilcup, Patel 
and Sharpe (GKPS) , which laid a foundation 
for a lattice QCD determination of this quantity 
with the Kogut-Susskind quark action. A num- 
ber of problems, however, are posed to warrant 
viability of the lattice result in the continuum. 
While some aspects have been examined previ- 
ously |^,|^, further studies are needed to estab- 
lish the continuum value of Bk- The statistics 
of GKPS are only 15 (/3=6.0) to 24 (/3=6.4) con- 
figurations: as a consequence the error of Bk at 
P = 6.4 is so large that an 0{a'^) dependence of 
Bk on the lattice spacing inferred from a theo- 
retical analysis |3| can not be ascertained. Fur- 
thermore the lattice size of 32"^ x 48 employed at 
(3 — 6.2 and 6.4 may be too small in the tempo- 
ral direction. We have therefore decided to carry 
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out a new set of simulations with improved statis- 
tics, incorporating advance in methods of analy- 
ses. Here we present a status report of this pro- 
gram which is being carried out on VPP500/80 
newly installed at KEK. 

2. Simulation method and run parameters 

We basically employ the formalism developed 
by GKPS for calculating the B parameter. Fol- 
lowing improvements are included in our simula- 
tions, however. 

On each configuration we calculate 2^^ = 8 
quark propagators corresponding to wall sources 
with a unit source placed at one of the corners of 
each spatial cube. This allows us to construct a 
wall source which creates only the pseudoscalar 
meson in the Nambu-Goldstone channel 1^. In 
contrast the even-odd source employed in previ- 
ous studies ||l|-^ also creates vector mesons that 
contaminate signals. The sources are placed at 
the two edges of the lattice, and the Dirichlet 
boundary condition is employed in the temporal 
direction. Gauge configurations are fixed to the 
Landau gauge. 
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Table 1 

Run parameters 



/3 


5.85 


5.93 


6.0 


6.2 


L^T 


16^32 


20^40 


24^40 


32^64 


#conf. 


60 


50 


50 


40 


a-i 


1.34(3) 


1.58(4) 


1.88(5) 


2.62(9) 


iriqa 


0.01-0.04 0.01-0.04 


0.01-0.03 


0.005-0.02 


msa/2 


0.0197 


0.0156 


0.0125 


0.0089 


fit 


10-20 


12-26 


14-24 


20-42 



Four-quark operators written in terms of 
Kogut-Susskind quark fields are extended over a 
2^ hypercube. We employ both gauge-invariant 
and non-invariant operators Q|, which differ by an 
insertion of appropriate gauge link factors for the 
former. The one-loop result for the lattice renor- 
malization factors is used to convert lattice 
values of Bk to those in the continuum, renor- 
malized with the naive dimensional regulariza- 
tion(NDR) scheme. 

A source of uncertainty in the conversion is the 
choice of the gauge coupling constant used for 
the renormalization factors. We take the mean- 
field improved MS coupling constant at the cut- 
off scale /i = 7r/a[ p[. Correspondingly we im- 
prove the operators [|10| through the replacement 
X a/moX and f/^ u^^Ufj,, where x is the 
Kogut-Susskind quark field, the link variable 
and Mo = P^/^ with P the average plaquette. 
Modifications of renormalization factors due to 
this replacement are also taken into account. 

We note that we have ignored the mixing 
of non-Nambu-Goldstone four-quark operators in 
the present analysis. Studies of this mixing ef- 
fect, which is believed to vanish in the continuum 
limit, are left for future work. 

Parameters of runs we have so far carried out 
are listed in Table |l|. Gauge configurations are 
generated with the 5-hit heatbath algorithm at 
2000(/3 < 6.0) or 5000(/3 = 6.2) sweep inter- 
vals. The physical scale of lattice spacing is de- 
termined from the p meson mass in the VT chan- 
nel. Values of /3 are chosen so that the lattice 
spacing is roughly equally spaced, and the lat- 
tice size is increased for weaker coupling to keep 
the physical lattice size approximately constant 
at La K, 2.4 — 2.5 fm. A similar consideration is 
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Figure 1. Comparison of lattice value of Bk for 
gauge non-invariant operator at {3 — 6.0 and 6.2. 
Data for j3 = 6.2 are shifted downward by 0.05 to 
avoid overlap. 

made for the choice of the fitting interval in ex- 
tracting Bk from kaon Green's functions. Quark 
mass is taken in steps of 0.01 for (3 < 6.0 and 
0.005 at /3 = 6.2. While we have evaluated Bk for 
non-degenerate strange and down quark masses, 
we shall restrict ourselves to the degenerate case 
in the present analysis. The strange quark mass 
yielding the observed ratio rriK/mp = 0.648 for 
this case is also given in Table |l|. We estimate er- 
rors by a single elimination jackknife procedure. 

3. Comparison with previous work 

In Fig. 1^ lattice values of Bk for the gauge non- 
invariant operator at /3 = 6.0 and 6.2 are plotted. 
At (3 = 6.2 our results are consistent with those 
of GKPS[^|j|] at the level of one standard devia- 
tion. On the other hand, our values at /? = 6.0 
are systematically lower than those of GKPS and 
Ishizuka et (IFMOSU) ||] obtained on the same 
size of lattice. 

In order to understand the discrepancy, we note 
that GKPS and IFMOSU respectively used 15 
and 10 configurations as compared to our 50, 
and that both groups employed the even-odd wall 
source, while we used the 2^ source. For a com- 
parison under the same condition, we divided our 
ensemble into 5 sets of 10 configurations and re- 
calculated Bk for each set. In addition we made 
a new calculation for 3 sets of 10 configurations 
using the even-odd wall source. The eight values 
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Figure 2. Lattice value of Bk at f3 — 6.0 calcu- 
lated for 10 configurations with 2^ (circles) and 
even-odd (diamonds) source as compared with 
previous results obtained with even-odd source 
(triangles). Solid lines are values for the full en- 
semble with errors indicated by dotted lines. 

of Bk thus obtained are plotted in Fig. || and 
are compared with those of IFMOSU and GKPS. 
Solid lines represent values obtained for the full 
ensemble for each type of source. 

We note that, while results for the full ensemble 
are consistent between the two types of sources, 
both errors and fluctuations of central values are 
significantly reduced with the 2^ source, demon- 
strating the advantage of this method. Since val- 
ues of GKPS and IFMOSU are consistent with 
ours for the even-odd source within one to two 
standard deviations, we conclude that previous 
results have overestimated the value of Bk at 
P — 6.0 due to insufficient statistics. 

4. Quark mass dependence 

Extraction of the physical value of Bk requires 
an interpolation of results to the quark mass cor- 
responding to the physical kaon. We find a sug- 
gestive result that the quark mass dependence of 
our data is well described by the form predicted 
by chiral perturbation theory, which we employ 
for the above purpose. For the case of degenerate 
strange and down quark treated in the present 
analysis, the formula is given by|pd[| 

Bk = B[1 - 3c{mqa) log(mga) -I- b{mqa)] , (1) 

where c — a/ {An f^^)"^ with a = (mKCi)'^ /mqa and 
b is an unknown constant. A representative ex- 
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Figure 3. Mean-field improved Bk(NDR, 2GeV) 
as a function of rriqa at (3 — 6.2 on a 32^ x 64 
lattice. Solid lines are fits to (Q) and dashed line 
corresponds to physical value of mK/mp. 

ample of the fit taking B, c and b as parameters 
is plotted in Fig. ^ together with the quark mass 
corresponding to mK/mp = 0.648 (dashed line). 
Several uncertainties exist, however, to take the 
agreement as evidence for the presence of chiral 
logarithm in our data. The pion decay constant 
estimated from the fitted values of c and the 
spectrum results for the slope a turned out to 
have a large value of 150-300 MeV with 10-20% 
error as compared to the experimental value of 93 
MeV. Also a quadratic form without logarithm 
fits our data equally well. 

5. Results for Bk 

We plot our results for Bk renormalized at 
/i = 2 GeV in the NDR scheme as a function of 
TTipfl in Fig. ^ both for gauge invariant (filled cir- 
cles) and non-invariant (open circles) operators. 
We observe that the data do not exhibit the the- 
oretically predicted O(a^) dependence[||, and are 
apparently more consistent with a linear behav- 
ior, at least over the range nipa = 0.57 — 0.29 (or 
1/a — 1.3 — 2.6 GeV from rrip) corresponding to 
/3 = 5.85 — 6.2. If we assume a linear form, we 
obtain ^^(a) = 0.4825(83) + 0.489(17)mpa and 
0.5097(83) -H0.464(17)mpa for the gauge invariant 
and non-invariant operators. 

We have varied the analysis procedure to see to 
what extent the values of Bk are modified: (i) In- 
stead of niK I JTip we examined the choice /rup 
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Figure 4. Mean-field improved ^/^(NDR, 2GeV) 
as a function of mpO, together with hnear fits 
(solid lines). 

for estimating the strange quark mass. This only 
resulted in a common downward shift of the four 
values of Bk by about 0.02, and hence does not 
change a linear dependence on the lattice spacing, 
(ii) We also replaced the mean-field improved MS 
coupling constant gjjg{-K/a) in the renormaliza- 
tion factors by the running coupling constant at 
jjL = 2 GeV evaluated by the two-loop renormal- 
ization group evolution starting from ctjjg{'!T / a) . 
This increases the value of Bx by about 0.03-0.05 
without, however, changing the linear behavior. 

Let us comment on the fact that our results 
for Bk obtained with the two types of opera- 
tors do not agree with each other and that the 
discrepancy does not become smaller toward the 
continuum limit. Possible sources of the discrep- 
ancy are an operator dependence of finite lattice 
spacing effects, which disappears in the contin- 
uum limit, and 0(g*) corrections in the renor- 
malization factors which are not included in our 
analysis. Naively one would expect the magni- 
tude of 0{g^) corrections to also diminish toward 
the continuum limit. This is not the case, how- 
ever, for Bk renormalized at a fixed scale /z = 2 
GeV, since the gauge coupling constant g(/i) gov- 
erning these corrections is independent of the lat- 
tice spacing. We have checked that increasing the 
scale n leads to a smaller discrepancy. Thus the 
difference of Bk for the two types of operators 
extrapolated to the continuum limit provides an 
estimate of the magnitude of 0{g'^) corrections. 



Finally we remark that an 0{a?) dependence 
is not completely ruled out from our data if there 
exists an appropriate 0{a?) or O(a^) term which 
just brings the curve to behave as if it is linear 
over a range of lattice spacing including that ex- 
plored in our simulation. If this rather excep- 
tional case happens, the value of Bk extrapolated 
to the continuum limit a = would increase by 
about 0.05. 

6. Summary 

In this report we have described our effort to- 
ward improving previous work on the kaon B pa- 
rameter with the Kogut-susskind quark action in 
quenched QCD. The results obtained at four val- 
ues of (3 show that Bk falls linearly with the lat- 
tice spacing over the range 1/a = 1.3 — 2.6 GeV. 
If the theoretical prediction of an 0{a?) depen- 
dence is valid, this means that such a behavior 
will become apparent only for smaller values of 
lattice spacing. We plan to carry out a simula- 
tion at /3 = 6.4 — 6.5 employing the lattice size of 
40^ X 96. 
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